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Abstract 

We show that the first non-linear correction to the linearised Poisson- 
Boltzman n (or DLVO) theory of effective pair interactions between charge- 
stabilised, co lloidal particles near a charged wall leads to an attractive com- 
ponent of entro pic origin. The position and depth of the potential compare 

favourably with rec ent experimental measurements 
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Recent direct measurements confirm the overall validity of the classic Derjagui n-Landau- 
Verwey-Overbeek (DLVO) effective pair potential between charge-stabili zed colloidal parti- 
cles in the bulk of an aqueous dispersion |T|], but give convincing evidence for an attractive 
component to the interaction when two co lloidal particles are confined in a slit |||3| or close 
to a planar wall cite4. In this letter we use a density functional analysis to show t hat 
such an attraction can be due to non-linear corrections to linearised Poiss on-Boltzmann 
(PB) theory which forms the basis of the DLVO potential. 

Mesoscopic colloidal particles carrying ionizable radicals on their surface build up elec- 
tric double -layers of microscopic co and counterions when dispersed in water. The effective 
forces between these "dressed" particles or polyions derive from a sum of direct and induced 
interactions. The direct interactions include excluded volume and Coulomb repulsions be- 
tween the bare surface charges, and van der Waals attractions; the latter are generally quite 
negligible, except at very high salt concentrations corresponding to the strong screening 
regime |(|. The indirect interactions are induced by the locally inhomogeneous distribu- 
tion of microions, which screen the bare Coulomb repulsion; these induced interactions are 
identified with the free energy of the microions in the "external" field due to any given con- 
figuration of polyions, and hence have electrostatic and entropic contributions. In practice 
the free energy of the microions, which depends parametrically on the positions {Ri} of 
the N polyions (assumed here to be spherical) is calculated within density functional the- 
ory (DFT) by minimising an appropiate free energy functional F with respect to the local 
densities p Q (r) of the co and counterions. To be specific, consider polyions of radius R, car- 
rying a total positive charge Ze {Z > 0), and monovalent coions (local density, p+(r)) and 
counterions (p_(r)). Z must be regarded as an effective charge much smaller than the bare 
structual charge. The "renormalized" charge accounts for the shortcomings of linearised 
PB theory (as used later in this paper) at very short distances from the polyion surfaces 
jTJ. Correlations between microions will be neglected, which amounts to working within the 
framework of mean- field, or Poisson-Boltzmann theory. The solvent (water) is merely treated 
as a continium of dielectric constant e ("primitive" model). Under these conditions the free 
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energy functional reduces to the sum of ideal and electrostatic contributions: 

F[p+,p_] = F td [ P+ ,p_] +F el [p + ,p_] (1) 
Fu\p+,P-] = E k BT [ Pa(r) [log(A^ a (r)) - l] dr (2) 

Fei = ~ / p c (r)*(r)dr (3) 

Z Jp 

where the A a are irrelevant length scales, p c (r) is the total charge density: 

N 

Pe(r)=EPex t W+e[p + (r)-p_(r)] (4) 

i=l 

Pelt( r ) is the "external" charge density associated with polyion i, p is the domain to which 
the dispersion is confined, and \I/(r) is the local electrostatic potential, which is related 
to p c ( r ) by Poisson's equation; the latter must be solved subject to appropiate boundary 
conditions on the surface bounding p. 

If the p Q (r) are slowly varying, the ideal contribution may be expanded in a functional 
Taylor expansion around the mean (macroscopic) densities n a . If the expansion is truncated 
after second order, which amounts to linearisation of PB theory, the optimum density profiles 
(which minimise the free energy functional), and the resulting electrostatic potential, reduce 
to sums of terms ("orbitals") centred on each of the N polyions 0: 

N 

Ap Q (r) = p a (r) -n a = £p«(r - R ; ) (5) 

i=l 

where p^(r) is constrained to vansh for | r |< R (excluded volume). Apart from a structure- 
independant "volume" term || the total potential energy of the "dressed" polyions reduces 
to a sum over all pairs of an effective pair potential: 

v eff (R b Rj) = e / pgt (r)*®(r)dr (6) 

Jp 

where ^^(r) is the electrostatic potential due to polyion j and its associated "orbital" 
p+\r) — p_ (r). In the bulk, i.e. far from any surface bounding p, v e ff reduces to the familiar 
DLVO potential 0; 
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, e// (Rii) = Z c 2 ff e 2 eXp{ £ ^ Rij} (7) 

where k—1/Xd is the inverse Debye screening length, n = [4n(n + + n_)e 2 /ekBT] 1 ' 2 , and 
Z e ff > Z is an effective polyion charge which accounts for polyion-microion excluded volume, 
resulting in p^(r) = 0, | r |< R. 

We have recently shown how t> e yy(Ri, Rj) is modified when the two polyions i and j 
are close to a planar surface, which separates the dispersion from a medium (e.g. glass) 
of dielectric constant e' (see fig.l). Using an integral representaion of ^-^(r) due to 
Stillinger [11], which takes proper account of image charges, we were able to show that 
within linearised PB theory, the effective repulsion between two polyions situated at an 
altitude h above the surface is enhanced relative to the bulk result (^), and decays like 1/-Rfj 
at large (horizontal) separations Rij. In fact, a careful asymptotic analysis of Stillinger 's 
expression for *&v>(s, z) (where s is the radial coordinate parallel to the plane, and z the 
vertical coordinate) shows that: 



-KT g 



KT ' o~ k(2H-\-z) 

+ —r + —2 — + °^ 



(8) 



where all distances (s,z,r = (s 2 + z 2 ) 1 / 2 and r' = (s 2 + (z + 2h) 2 ) l l 2 ) are relative to the centre 
of the polyion j, (or its image j') and e r = e'/e. The resulting v e jf easily follows from eq. (P) 
and has the same asymptotic functional form. 

It is thus clear that linearised PB theory cannot account for the attractive component of 
the effective pair potential between polyions, which is clearly apparent in recent experimental 
determinations of f e //(Ri,Rj) in confined geometries [2-4] (for a possible interpretation of 
the experimental data, see however [Q). The conclusion from linearised PB theory is 
independant of any surface charge carried by the planar boundary. To make progress, the 
full (non-linear) PB theory should be used. This can only be done numerically and represents 
a formidable task, in view of the abscence of simplifying symmetries. We hence resort to a 
perturbative approach, whereby the expansion of the ideal entropic terms in eq.(^) around 
the uniform densities n a is taken to third order, yielding: 



Ap a (r) 



F id [p + ,p„] = F td (n + ,n„) + k B T^2 I APa(r) [log^n,,) + -^Idr (9) 



2n a 



-^E/^r + 0(A^) 



linearised PB theory amounts to retaining only terms of up to quadratic order in Ap a , and 
leads back to the linear superposition @ of "orbitals" of the form given in and JT1] . In 
a density functional perturbation theory similar in spirit to that used by in references || 
and Jl3| to determine triplet interactions in the bulk, this linear superposition of orbitals is 
then substituted into the second integral, involving cubic terms in Ap a (r), in the functional 
expansion (g). 

Since we are interested in the effect of a charged planar surface (carrying a uniform 
surface charge a) on the effective pair potential between two polyions (say i = 1 and 2), the 
relevant superposition of "orbitals" , to be substituted into the perturbation term is; 

A Pa (r) = pP> + p<J> + ;<* = +,- (10) 

where p£'(r) is the planar "orbital", associated with the charged surface, as calculated 
within linearised PB theory, i.e. 

p<g){z) =oKe-< z ^ (11) 

Introducing the local charge density in p: 

Ap(r)=Ap + (r)-Ap_(r) (12) 
= E[p?W-P (i) W]=EP W (r) 

i i 

the lowest-order non-linear correction to the effective pair-potential between polyions 1 and 
2 reduces to: 

Af e //(Ri,R 2 ) = k f ~ / p«(r)p( 2 )(r)^)(r)dr (13) 



k f {n Z - " 3 +) / P«(r)p (2) (r-R 12 )p M (r)dr 

n_ + n + V Jp 



5 



where p^(r) is taken from the linearised theory 0,11] and 



calculated from eq.(|Tl|). The correction in eq.(13) is the lowest order (cubic) term which 
takes into account the presence of a charged surface in the vicinity of the two interacting 
polyions 1 and 2. There is good reason to believe this correction should be well behaved since 
for reasonable separations and distances from the wall no two terms will simultaneously be 
large at any one position in space. Linear theory (which includes terms only up to quadratic 
order in Fid), on the other hand, accounts for the dielectric discontinuity at the surface, which 
reflects itself in the introduction of image charges |]lO , [lT| , but is "blind" to the presence of 



an electric double layer (i.e. an inhomogeneous distribution of co and counterions) in the 
vicinity of the charged surface. An explicit albeit lengthy, expression for Au e // in k-space is 
obtained by taking the Fourier transform of eq. ([13|) and invoking the convolution theorem. 

The inverse Fourier transform, to revert to r-space, can be reduced to a ID quadrature 
which must be carried out numerically and leads to a result of the form; 

Av eff {h u h 2 , R 12 ) = -k B T Zeff<JK (n ~~^ )e -f(K h 2 , R 12 ) (14) 

(71+ + Tl—) 

where f is a dimensionless function of the arguments hi,h 2 , the heights of the polyions, and 
i?i2, the distance between them. In the case where (hi — h 2 — h), f is given by, 
R v_ r Uo(lKR 12 )dl { l + g(l)e-^ 1/2 « h 1 - eg-g+g^ 

/(/ ^ i2) -y p+i [ 2(/2+ i)V2 + i + 2(p+ 1)1/2-1 (i5) 

+ 2 (J* - l) e-^^g{l) + ^^^ (1 + 9(l)e Kh )"j 



where, 



(12 _j_ 1N1/2 _ 7 

9(1) 



+ 1)1/2 + £p J 

Note that the charge <r may be of the same sign as that of the polyions, or of opposite sign. 
In the latter case the sign of (n_ — n+) will depend on the system. In the canonical ensemble, 
where the number of ions is well defined by charge neutrality, (n_ — n + ) will be negative 
for a sufficiently high surface charge; while in a semi-grand canonical scheme, where the 
system is in contact with a reservoir of co and counterions, (n_ — n + ) will be determined by 
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the chemical potential of the co and counterions in the reservoir. The asympototic form of 
Av e jf for large Rij is, 

Zt ff g K *(n--n + )e-* h ( 4e-* R « 2e^ h (l - e^)e^_ 
Av eff - -k B T (n+ + n _ )3 I —3— + ^ ( 16 ) 



2e r e~ Kh 



(7 - Ae' Kh ) ^ 



We have calculated t> e // for two polyions situated at an altitude h above the charged plane, as 
a function of the horizontal distance d = \ R42 |. Examples of the total effective pair potential 
(including Av e ff) are plotted as a function of d in Figures 2 and 3. in Fig. 2 the calculated 
potentials are shown for several altitudes h, at a fixed surface charge a = 16000e//im 2 .This 
surface charge a must once more be considered as an effective charge; the tightly bound 
counterions of the "Stern layer" will not contribute to Av e ff m eq.(13). This surface charge 
is limited to a maximum value of 21000e//im -2 by the slit width and the screening length 
Xd quoted in ref. [[|; the contact-potential corresponding to o = 16000e//im~ 2 is 3 in 
dimensionless units within linear theory . The physical conditions (polyion size R, charge 
Ze and concentration n; monovalent salt concentration n s in a closed system) are as close as 
possible to the experimental conditions of ref. Q (this reference does not specify the surface 
charge a of the glass wall which is difficult to measure). As the polyions move closer to 
the surface, an attractive potential well deepens rapidly and shifts to smaller separations 
d. Both the depth and position of the potential minimium agree qualitatively with the 
experimental data ||, although the calculated attraction appears to be of shorter range. 
These features vary rapidly with altitude h; at the lower altitudes, when (h — R) becomes 
comparable to Ad (which is of order of 0.28 fim under the present conditions), higher-order 
terms in the expansion of Fid are expected to become important. 

In Fig. 3 the same total effective pair potential is plotted versus d for several surface 
charge densities a. The variation with a is seen to be again very rapid; as a increases, the 
well depth increases and its position shifts to smaller values of d. Note that if a is of opposite 
sign to Z\ Av e ff may still be attractive, since (n_ — n + ) may be negative. 
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In summary we have shown shown that while linearised PB theory invariably predicts 
a purely repulsive effective interaction between polyions in the bulk or near a surface, per- 
turbative inclusion of the lowest order non-linear term into the ideal contribution to the 
free energy functional leads to an attractive component of the effective pair potential. The 
depth of the resulting potential well is of order fcgT under physical conditions comparable to 
those of recent experiments [2-4]. Although the attractive component vanishes with a, the 
effect is of partially entropic origin, and is due to a competition of electrostatic and entropic 
couplings between the "orbitals". Since this paper was submitted, a numerical solution of 
the full non-linear PB theory in a different geometry (two colloidal particles inside a charged 
cylinder) has been published |TJj . This numerical work also predicts an attractive well in the 
effective polyion-polyion pair potential, in qualitative agreement with the present analytic 
study. 
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FIGURES 
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r'(s,z-2h) 



FIG. 1. Geometry of particles near a planar surface. The suspension, of dielectric constant e, 
is confined to z > —h, while a material of dielectric constant e' occupies the half-space z < —h and 
carries a surface charge a. The dotted circles represent the image charges. 
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FIG. 2. Dimensionless total (Linear + Perturbative solution) effective pair potential, 
v e ff(h; dj/ksT versus distance d between polyions situated at the heights shown in the legend, mea- 
sured in /im. The polyions have a bare valence Z=5000, and a diameter D = 0.65/um. T = 298-ftT; 
e = 78 and e' = 6.3; n s = 3 x 10~ 7 M and the wall carries a charge of 16000 e/jim 2 
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FIG. 3. Same as in figure 2 but this time the height of the particles is fixed at 1.5 fim and 
the surface charge is varied as shown and is measured in e//im 2 ; the Debye length is fixed at 
Xd = 0.28/xm, so that the salt concentration n s varies with a 
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